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The Dvali-Gabadadze-Porrati model introduces a parameter, the cross-over scale r c , setting the 
scale where higher dimensional effects are important. In order to agree with observations and to 
explain the current acceleration of the Universe, r c must be of the order of the present Hubble radius. 
We discuss a mechanism to generate a large r c , assuming that it is determined by a dynamical field 
and exploiting the quantum effects of the graviton. For simplicity, we consider a scalar field ^ with 
a kinetic term on the brane instead of the full metric perturbations. We compute the Green function 
and the 1-loop expectation value of the stress tensor (T M „) of *P on the background defined by a 
flat bulk and an inflating brane (self-accelerated or not). We also include the flat brane limit. The 
quantum fluctuations of the bulk field ^ provide an effective potential for r c . For a flat brane, the 
1-loop effective potential V c tt(r c ) is of the Coleman- Weinberg form, and admits a minimum for large 
r c without fine tuning. When we take into account the brane curvature, a sizeable contribution 
at the classical level changes this picture. In this case, the potential can develop a (minimum) 
^ ' maximum for the (non-) self-accelerated branch. 
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I. INTRODUCTION 



. One of the most startling observations in cosmology in the last few years is that our universe is currently in a 
f*** ' phase of accelerated expansion P, . One possibility to explain this phenomenon is to modify the gravitational 
sector at large distances. An interesting realization of this idea was proposed by Dvali-Gabadadze-Porrati (DGP) 0, 
and arises within the Brane World scenario. In the DGP model, there are two curvature terms, one in the bulk and 
the other on the brane. The ratio between the respective Newton constants 1/M 3 and \/m 2 P defines an important 
quantity, the crossover scale r c . At distances shorter than r c , gravity behaves as in four dimensions whereas for larger 
distances it behaves as in five dimensions. 



This model admits a 'self-accelerated' solution, where the brane accelerates at late times even with vanishing brane 
tension or bulk cosmological constant Q. In this solution, the acceleration arises as a purely gravitational effect, 
and the Hubble rate is On the other hand, r c must be of the order of the present Hubble radius in order not 

<— i ' to conflict observations 

IE HHH (see also 0] for a recent review on the phenomenology of the DGP model). Thus, 
this model provides an appealing explanation of the current acceleration of the Universe. There exists an abundant 
literature concerning the stability of the self-accelerated solution |lfj. III! Il2l Il3l Il4l Il5|. So far, there seems to be 
agreement that perturbatively, the DGP model has ghosts for the self accelerating solution. However, this mode is 
strongly coupled, so it is not clear if it remains a ghost at the full non-perturbativc level. Here, we shall leave this 
issue aside and, instead, we discuss another essential ingredient for the model to account for the current accelerated 
expansion, namely, whether it is possible to have a large r c naturally. 

One argument to obtain a large r c is that the quantum effects from the numerous matter fields on the brane would 
generate the localized kinetic term for the graviton with a large coefficient To this end, though, one would 

need an enormous number of fields. Another argument is that the cutoff on the brane and in the bulk could actually 
be different \vf\. We shall also mention that in the model of the corresponding cross-over scale is naturally large. 
Here, we shall take the approach that r c is determined by some (scalar) field, whose dynamics drives r c to a large 
value. Some mechanisms of this sort were already introduced in [lg. Our aim is to elaborate on this idea. We shall 
concentrate on the possibility that the quantum effects of the bulk fields (the graviton) provide such a mechanism, 
given that this requires a minimal extension of the model. The main assumption is that the Newton constant in bulk 
(or on the brane) is promoted to a Brans-Dicke-type field <&, which determines the value of the crossover scale r c . 
The idea is then that the quantum effects from the graviton induce an effective potential for $ that drives r c to large 



* pujolas@ccpp.nyu.edu 



2 



values. As we shall see, the 1-loop effective potential induced by the bulk graviton is of the form 

rri-loop, s _ OL-fi log(Ar e ) 
V cS I.M - 4 

c 

where j3 is a constant of order 10 -3 (see below for details), A the cutoff of the theory and a is an order-one constant. 
This potential has a minimum at r c ~ e"^/A, which can be naturally very large. This is a very appealing way to 
generate a large r c , but is not the end of story. The reason is that in the most interesting case when the brane is 
accelerating, there is a large contribution at the classical level due to the curvature terms in the bulk and on the brane. 
In the simplest situation, the minimum of the potential occurs at a value of r c too small for phenomenological purposes. 

As mentioned before, here we discuss the quantum fluctuations in the DGP model. Specifically, the aim is to 
estimate the vacuum fluctuations of the graviton on the background defined by a 4D inflating brane embedded in a 
flat 5D space. We shall consider the two possible branches described in £|, and also the case when the brane is flat. 
To simplify the analysis, instead of the full graviton we shall consider a 5D massless minimally coupled scalar field ^. 
This replacement is motivated because in an appropriate gauge, the spin-2 metric fluctuations obey the same equation 
as a massless field, so the analysis of the Green function do not differ much. In practice, then, we are ignoring the 
brane bending mode of the true DGP model. The discussion of this sector is quite more complicated because it is 
strongly coupled 0, [n| , aside from the fact that it is a ghost in the self-accelerated branch [ToL [Til IT^ . Il3| . On the 
other hand, quantities like the vacuum expectation value of the stress tensor (T M1/ ) may certainly differ for the true 
graviton and the scalar considered here, but still we shall restrict to the latter for illustrative purposes. 

On the more technical side, we shall mention that the above result for the potential seems badly divergent in 
the r c — > limit. As we show in the body of the paper, this is only an artifact of the renormalization procedure 
used. Schemes based on analytic continuation (dimensional regularization, £— function) lead to the expression above. 
Using a more physical scheme, for instance introducing a cutoff, one unveils divergences of the form A k /r^ k with 
k = 1,2,3,4, which are implicitly 'renormalized away' by the schemes based on analytic continuation. Including 
these terms, the r c — ► limit is regular as expected. Still, for large r c , it makes sense to use the usual rules of 
renormalization in terms of the effective field <p = l/r c . As we shall see, this is already suggested by the fact that the 
classical ('boundary') effective action [T^j for the DGP model is analytic in l/r c (for large r c ). 

The quantum effects from bulk fields with brane-localized kinetic terms have been considered previously in a 
slightly different context. The authors of computed the Casimir energy present in a system of two flat branes, 
and showed that when the weight of the localized kinetic terms are different, the Casimir force could stabilize 
the distance between the branes L. The use of the Casimir effect to stabilize the size of the extra dimension 
was first discussed in (53, El, and in the brane world setup in [2^, l23l l24j . The work of is relevant here 
because in the one brane limit (when one of the brane is sent to infinity), their result must be compatible with the 
effective potential found here in the limit when the brane is flat. The L — > oo limit of the result quoted in |l9| 
vanishes, apparently suggesting that the potential V e g(r c ) is trivial, which is of course not the case. The reason for 
this apparent discrepancy is that reference 0] deals with the L— dependent part of the potential (the Casimir energy). 

This paper is organized as follows. In Section [H] we introduce the scalar version of DGP model. In Section ITTT1 
we work out the fluctuations of the model in the background defined by the de Sitter brane (in both branches). We 
compute the Green function in Sec IIII Al and (T M „) in Sec IIII Bl In Sec IIII CI we discuss the limit when the brane 
is flat, and Sec IIII Dl deals with the stress tensor on the brane. The reader not interested in the technical details 
might skip these sections, and go to Section llVl where the effective potential for r c is derived, and Sec II V Al where 
the mechanism to obtain a large r c is discussed. 



II. THE DGP MODEL 



The action for the DGP model is 



S = J d b x^fg~ M 3 R {5) + J d A x\fh(m 2 P R (i) -T) , (2.1) 

where M and mp are the 5 and 4-dimensional Planck masses, R(S) is the bulk Ricci scalar, g and h denote the 
determinants of the induced metrics in the bulk and on the brane, R(a) is the Ricci scalar of the brane geometry, and 
r is the brane tension. This model admits two solutions with a line element of the form £| 

ds 2 = dy 2 +a 2 {y)ds 2 (i) (2.2) 
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where ds 2 ^ is the metric on a de Sitter space of unit curvature radius and with 

a(y) = H- 1 ± \y\ (2.3) 

where H is the Hubble constant on the brane (located at y = 0). Throughout this paper, the upper sign corresponds 
to the self accelerated branch (or simply the '+' branch) Q. The — or 'normal' branch can be visualized as the 
interior of the space bounded by a hyperboloid (the brane), whereas the + branch corresponds to the exterior. For 
each branch, the Israel junction conditions require that 



±H = r c H 2 , (2.4) 

12M 3 ' V ; 



where we introduced the cross-over scale 



~ 2Af 3 ' 

Eq. H2.4|) fixes H in terms of r c and r. In the literature, the term 'self-accelerated' is usually reserved for the solution 
with t — (for which Hr c = 1). Here, we shall consider all values of Hr c , which can be obtained by appropriately 
choosing r and we call 'self accelerated' any solution in the + branch. 



A. A dynamical cross-over scale 



constant on the brane or in the bulk so that ($5) = M 3 , and ($4) = m%. In terms of them the effective cross-over 

(2.6) 



We now extend the original induced gravity term to the most general Brans-Dicke-type modification of the DGP 
model, 

S = J d 5 Xyfg ($ 5 R {5) + kin) + / d^xjh ($ 4 %) + kin) (2.5) 

where $5,4 are unrelated bulk and brane fields and we have omitted their kinetic terms, as we are interested in the 
configurations with constant fields only. The important point is that each BD fields control the value of the Newton 

? 4 ) = nip. 

scale is give by 

r =-?! 

Tc 2$ 5 

The equations of motion in the model l|2.5[) become 

V eS = -12<J> 4 iJ 2 ± 12$> 5 H + t + 5V(<f> 5 , $ 4 ) = , (2.7) 

together with V^' ff = 0, where a prime denotes here differentiation with respect to $5 or <5> 4 . In Eq. (|2.7|l . 5V 
represents the potential responsible for fixing them (which might include bulk and brane contributions). Note that 
the condition l|2.7|l is nothing but the extension of the Friedmann equation l|2.4|l . The first two terms in l|2.7|) are 
classical contributions due the brane and the bulk curvature terms in the action. The former induces a negative 
runaway potential for $4, while the latter is positive for the self-accelerated branch. As we shall see, for large r c , the 
1-loop effective potential is analytic in l/r c cx 1/<I>4. Thus, it seems difficult that the quantum effects can induce a 
minimum at a large $4. From now on we shall consider that $4 = rap is fixed, and instead we try to 'stabilize' $5 
using the quantum effect from the graviton. Note also that this choice automatically leads us to the four dimensional 
Einstein frame, which is the physical one. See |3(j for a recent discussion of the cosmology of a DGP scenario with a 
dynamical $4 (and constant $5). 



B. Scalar toy model 

We are interested in the fluctuations of the graviton propagating in the background 12. 2|) . One can show that the 
spin-2 polarizations of the graviton perturbations obey the same equation of motion as a massless minimally coupled 
scalar field (see e.g. |12|) with a kinetic term on the brane. Hence, we shall simplify the problem and rather consider 
a scalar field (even under Zi symmetry) propagating on the background l|2.2|) . Its action is 

S =-\f d 5 x^g(d^) 2 - J d 4 xVhr c (d^ f , (2.8) 

where, g and h denote the determinants of the background metrics on the bulk and on the brane (|2.2|) 
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C. Classical boundary effective action 



Let us now discuss how the cross-over scale appears in the effective four dimensional theory. Consider for simplicity 
the solution with a flat brane, and we have in mind a situation with large r c (compared to any other scale in the 
problem). A simple way to obtain the effective theory of the model (|2.8|1 is to compute the effective boundary action. 
Following [l0|. once some data on the brane are specified for the bulk field ^lo = ip, the equation of motion in the bulk 
is solved by 4" = e~ yA ?p, where A = (— Or^) 1 ' 2 and Ou\ is the 4D D'Alcmbertian operator. Inserting this solution 
back into the bulk action, we can integrate it explicitly and obtain a surface term. Thus, the boundary effective action 
is 



S eS = J d A xVh $ (r c D (4 ) - A) 



For large values of r c , the four dimensional term dominates. Hence, it is preferable to express this as the action for 
a 4D canonical field. Restricting our attention to constant configurations of r c , we can always rescale the field as 
tp = \ / 2r c if; in the canonical form, so the effective action is 

S ee =^[d A xVh(yn (4) $-v$Aty, (2.9) 

where we introduced 

1 

tp = — . 

r c 

This effective action describes the resonant or metastable mode of the DGP model. This mode propagates as in four 
dimensions during a typical distance r c after which it behaves as in five dimensions. Put another way, the decay width 
into five dimensional modes is is of order l/r c . Hence, at distances shorter than r c this 4D effective representation is 
accurate. 

The point that we want to emphasize here is that in this regime, the effective action (|2.9|1 is analytic with respect 
to ip. This is a rather special feature of this model. As we shall see, all the explicit expressions computed below are 
regular in the tp —> limit. Finally, let us note that in terms of the model of Sec III Al we have tp = 2<I>5 /mp, so we 
can interpret ip as the bulk Brans-Dicke field (restricted on the brane). 

Before switching to the explicit computation of the Green function, we shall briefly discuss the generalization of 
the boundary effective action to the curved brane case. It is straightforward to show from Eqs. (|3.1[) and 1|3.3|) 
that when the brane is inflating, the solution for '5 in the bulk with data on the brane specified by ip is given by 

= a~ 3 / 2 e - ' z 'PV ; : where z is the conformal coordinate (dz — dy/a), p = (— D( 4 ) + (3/2) 2 ) 1 / 2 and now D( 4 ) is the 
covariant DAlembertian in de Sitter space. Using the same logic as above, the boundary effective action in this case 




r eff = / d*xVh i/> r c D (4) =F -H — a I — - - DM) V ■ (2-10) 



From the analytic structure of this action as a function of D(4), we can already identify the main features of the 
effective theory. For both branches, there is a branch cut starting at D( 4 ) = (3H/2) 2 , which corresponds to the 
continuum of Kaluza-Klein modes. Aside from this, we can also see the presence of one lighter mode. Performing 
the same rescaling of ip as above, we can identify the mass of this mode as the non-derivative term in l|2.10|l . For the 
non-self-accelerated branch, this mode is massless. As we will see in short, this agrees perfectly with the Kaluza-Klein 
decomposition. For the self-accelerated branch, we read from (|2.10|l a mass-squared given by 3H/r c , which agrees 
with the analysis of the KK spectrum l|3.13fl to leading order (in the large Hr c limit). 

III. QUANTUM FLUCTUATIONS 

In this Section, we shall take the bulk spacetime to be (n + 2) dimensional. The Klein-Gordon equation on the 
space ()2.2|) . 

[D (n+2) + Hr c S(z)D (n+1) ] V = (3.1) 
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is separable into modes of the form 

H>(z,x»)=J2u p (zWp(^) > (3-2) 
p 

In the previous equation, □( n+2 ) and Q( n +i) are the bulk and brane D'Alembertians. For clarity, we omit the set of 
indices labelling the angular momentum of the four dimensional modes ip p (x 11 ). The (n+ 1)— dimensional modes obey 
[□( n +i) — m 2 ] tpp = and represent four dimensional fields with mass 

m 2 = [{n/2) 2 +p 2 ]H 2 . 

The mass spectrum is determined by the radial equation, which can be cast in the Schrodinger form 

[~d 2 + V(z)] U p = p 2 U p (3.3) 

with Up — a n l 2 Up, and 

n» = - t - ^ f> - H " 'tt ) 'W ■ <"» 

where a prime denote differentiation with respect to z, and we introduced the conformal coordinate dz = dy/a. 
Hence, the modes are combinations of e ±lp ' z ', and the spectrum is composed of a continuum of Kaluza-Klein modes 
with masses m KK > nH/2, plus possibly one discrete (or 'bound') state with mass m\, s < nH/2. 

In terms of the mode decomposition, the boundary condition on the brane reduce to 

[d z -v±(p)]U P \ 0+ = (3.5) 

where |o+ stands for lim 2 ^o with z > 0, and 

v ± (p)^±l-{p 2 +(l) 2 )Hr c . (3.6) 

Note that (|3.5|) is similar to the Dirichlet boundary condition for high KK mass modes, whereas it recovers the usual 
Neumann-type condition for light modes. This is indeed the key ingredient to obtain a four dimensional behaviour at 
short distances in the DGP model, since it implies that high energy modes are suppressed on the brane. 
We shall normalize the modes according to 

tf 2 ) = JdSfi (*td^ 2 - * 2 ^*D + 2 r c J da e i (*I0 P *2 ~ *2^**) lo (3.7) 

where is a future-directed timelike unit vector and <i£ {da) is the volume element of the surface S (a) normal 
to £ M in the bulk (on the brane). For any pair of modes satisfying the equation of motion (|3.1|l . the norm 13.7(1 is 
independent of the choice of E. We shall take £ M to be along the time coordinate of the observers on the brane. It 
is convenient to introduce the Rindler coorinates {a,i}, related to the rectangular coordinates in the bulk {T, X} 
through 

R = a cosh t 

T = asinh< (3.8) 

where R = |X|. In terms of these coordinates, £ M = Sp/a and the metric on the de Sitter slices y = const are 
^ s (n+i) = a [~dt 2 + cosh 2 t dfi 2 ] where dQ 2 n ^ is the metric on an n— dimensional sphere. 
For a mode of the form ty p j = U p ip Pt i where I is a collective index of the dS part, we have 



(*p,i,*p,i) = J dza n (z)(l + 2r c HS(z)) \U p \* <M ( „ +1) , (3.9) 

where (tpi, i/ , 2)(„ + i) is the usual four-dimensional Klein Gordon product, J da £ M i (ipldp,ip2 — ^d^tp*). The continuum 
of KK modes is normalized as ('J™, ^^/) = 5(p — p')6i,i', which leads to 

dz(l + 2Hr c 5{z))\U KK \ 2 =6(p-p') , 
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while for the bound state we shall impose 

f-OC 

dz(l + 2Hr c 5(z))\U bs 



A bound (or 'localized') state is a mode with pure imaginary p — iX and A > 0, so that its wave- function Ub s ~ e A ' 2: ' 
is normalizable. The boundary condition implies that A must satisfy 

-A-i/±(£A)=0 (3.10) 

where v± is given in (|3.6|1 . The two roots are, for each branch, 

= T 2 

T) I 

A 9 = ±- (3.11) 

2 Hr c y ' 

It is convenient to introduce the notation A> and A< as the larger and smaller of Ai^ . Only A> is positive, so this 
root corresponds to the bound state in both branches. Its wave-function is U\, s — N bs e~ x>z with (see l|3.7t ) 

_A> A; 

1 + 2Hr c A> nHr c =p 1 



N l = ; | orL x = rTTI^T ■ ( 3 - 12 ) 



In the — branch, the bound state is massless (A> = n/2). In the self accelerated branch, it is normalizable only for 
Hr c > 2/n, and its mass is 

2 nHr c — 1 

TO M+) = ~2 ' ( 3 - 13 ) 

' c 

The behaviour of this mode as a function of r c is as follows. For Hr c 2/n, the bound state is light. For smaller r c , 
it acquires a mass until at r c = 2/nH it degenerates with the KK. At this point, this becomes unnormalizable and it 
should be rather described as a resonance, or quasi-normal mode. Its mass acquires an imaginary part p = po — iT that 
is interpreted as the decay width into KK modes I n this model, this actually turns out to be a purely-decaying 

quasi-normal mode (po = 0), with decay width T = l/r c — 3H/2. 

The other mode, at A<, is not normalizable, so strictly speaking it is not included in the spectrum. However, one 
can think that it corresponds to a resonance, or equivalently that the whole KK continuum in some sense behaves 
like this mode. For the — branch, its decay width is l/r c + 3i//2, and for the + branch it is 3H/2. Note that in the 
latter case, we can have two different resonances for Hr c < 2/n. The spectra of normalizable and resonant modes is 
sketched in Fig. (1). 

A. The Green function 

Before starting the computation of the Green function on the backgrounds (|2.2|l . we shall mention that for a 
massless minimally coupled field, there is an ambiguity in the definition of the Green function in an overall additive 
constant. The shift symmetry ^ — > ^ + const allows us to 'gauge away' any constant contribution. This turns out 
to be important for the non-self accelerated branch because in this case there is a massless 4D mode, and as is well 
known massless minimally coupled scalar fields enjoy an (IR) divergent Green function in spaces with closed spacelike 
sections. In any case, and precisely thanks to the shift symmetry, this divergence does not appear in the stress tensor 
(Tpu). This issue has been discussed at length in four dimensions an d in the brane- world context j2f|. These 

explicit computations further show that both (Tj,„) and the Green function, are everywhere regular. 

We can compute the (Hadamard) Green function as a mode sum, 

G (l) = QKK + Qbs ^ with 

G b ° = fi(A>)W te (^( Z ')G<;» >(ds) and (3W) 
G KK = J o d P U™{z)U™{z')G% S) , 



1 To be precise, here we are assuming that Hr c is large enough 
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+ Branch 



Bound State X i3/2 - i/(H r c ) 
(or resonance) 



Resonance 



t 



KK continuum 
-i3/2 



(a) 



■ Branch 



Bound State X i 3/2 



Resonance 



KK continuum 



-i3/2 - i/(H r c ) 



(b) 



FIG. 1: We represent the spectrum of modes in the p plane, related to the mass by m 2 = (w 2 /4 + p 2 )H 2 . The continuum 
of KK modes have p real and are represented by the thick black line, both for the self-accelerated branch, (a), and for the 
non-self-accelerated branch, (b). The crosses denote the positions of the other two discrete modes in the spectrum for some 
representative value of Hr c . In the upper half plane, these modes are normalizable and correspond to the bound state. In the 
lower half plane, they are not normalizable, and represent a resonance. The shaded line indicates the possible values that the 
pole at p = i\2 can take as Hr c is varied from to oo. Note that in the self accelerated branch and for Hr c < 2/n there is no 
bound state. Rather, there are two resonances. 



where G^ dS ^ denotes the massive Bunch-Davies Green function in (n + 1)— dimensional de Sitter space. The 0(X > ) 

factor in G bs ensures that it is only included when it is normalizable. 

One easily finds [25| (for the 'right' half of the space, z > 0) that the normalized radial KK modes are 



w = y ^(zhZ+w^) ( cos + \ sin <p*>) ' ^ 

It is convenient to introduce the subtracted Green function 

G (1) = Cfi> - G2L . ( 3 - 16 ) 

where G^j ink is the Green function in Minkowski space, as it is free from ultraviolate divergences in the bulk. In 
terms of the mode expansion, 

G (1) = G bs + / dp [u; K (z)U™(z>)] sub G% s) (3.17) 

J — OO 

with 

[U« K (z)U™(z')] sub = £zifV<"+*> . (3.18) 

One can now perform the p integral in l|3.17|l closing the contour in the complex p plane by the upper half plane and 
summing the residues of the enclosed poles. Using 

P — = -l-2i±±±±(-± h-] (3.19) 

p + iv A> — A< \p — zA> p — i\ < J 

we see that the only pole from \U^ K (z)Up K (z r )\ sub contributing to (|3.17|) is at p = iA>, the bound state location. 
Given that 

A> +A< = 1 

A> — A< nHr c =p 1 ' 
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we find that this pole in the contribution from the KK modes equals the contribution from the bound state G bs (see 

■*(!) 

r p(dS)- 



Eq. H3.12|) ~). with the opposite sign 2 . The remaining poles in (|3.17|) arise from G^ 1 ),^. Performing the sum over 



them, (see the Appendix for details) we find 

G (1) = -Go - 2 ^ >+ ^< (A> G x> - A< G x< ) (3.20) 

where 

{ 2fj n / i \ 



G ° ~ nS (n+1) \1 + (aa'H 2 ) 2 -2aa'H 2 cos () '''' 2 ' ! 



G * ^ Hh^Vrh^- (-) 

Here, C is the invariant distance in dS space and <S'(„ +1 ) = 2n 1+n ^ 2 /r(l + n/2) is the volume of a unit n + 1 dimensional 
sphere, and 



oo oo 



F 1 (a,b 1 ,b 2 ,c;z 1 ,Z2)=}^}^ uaT, W ^ 3 - 23 ) 

fe=0 1=0 [ ,k+l 

is the Appell Hypergeometric function of two variables (here, (a)k = T(a + k)/T(a)). 

In the — branch A> = n/2, so it is clear from (13.22(1 and l|3.23[) that the A> term in (|3.20() induces a divergence, 
as previously anticipated. This divergence arises from the constant mode of the bound state, which is massless in 
this branch. It is well known that massless scalar fields in de Sitter space produce this type of infrared divergence. 
See |33l l34j for a discussion in dS space, and pl| for a similar brane model. Here, we will simply ignore all constant 
contributions to the Green function, given that they do not contribute to any observable associated with '5. This is 
a consequence of the shift symmetry of the model, \P — > + ccmst, which in other words allows us to gauge away the 
constant contributions. Hence, the leading order behaviour of H3.20[l near the light cone (aa' — > 0) is 

G W (y,y',x,x') ~ H n+2 aa' cos C . 



For a source point x' on the brane with t = 0, and x close to the horizon (U = R — T = 0), then £ is large and pure 
imaginary, so we have G ~ jjn+ij' ^ w } 1 i c ] :l j s perfectly regular on the horizon. 

For the + branch, we find that at infinity (aa' — > oo) the leading order terms in l|3.20|l cancel each other, and we 
have 

G (1) (y,y',x,x) ~ H n (aa'tf 2 )- ( ™ +1) cos ( . 

The disturbance generated by a source on the brane at t' = 0, at spatial infinity decays slightly faster than for a flat 
brane, G ^ ~ 1/ (Hy^ n+1 ^). See the Appendix for the form of G 1 in the bulk in the coincidence limit. 

B. Vacuum expectation value of the stress tensor 

The stress tensor for a classical field configuration can also be split into bulk and a surface parts as 



rpbrane n 

J-ij —tie 



1 

2 

6(y), (3.24) 



2 Note that any normalization of the modes other than 13.71 would not lead to this cancellation, with the corresponding wrong short 
distance behaviour of the Green function. 
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where h%j is the induced metric on the brane, and the equation of motion has been used. In point splitting regular- 
ization, the v.e.v. (T^) is computed as 3 



{T,u) bulk - ~ hm {d»dl-~g llv g x °d x d' < \G {1 \y,x ] y':x') and 

l y',x'^y,x { Z ) 

^brane = ^ ^ ]j m f Q.Q'. _ * h .. h kl d k d'\ (0, X; 0, x') , (3.25) 



where d' = d/dx 1 ^. 

For the — branch, we find 



, T yxbulk _ rrn+2 ) * _ X L_ ( 1 ( flg ) )" 1 fo 2 ft\ 

V y/ 2S (n+1) \ {\-{aH) 2 ) n + 1 l + nHr c (aH) 2 (1 - (aH) 2 ) n 1 ' ; 



nHr c + 2 „ 2 (n+l+l/Hr ) ^ / . , . 1 



{aH) B {aH) 2 \n + 1 + 



Hr c (l + nHr c ) y ' y ' \ Hr ( 

From the local conservation of the stress tensor, the remaining components (T 1 *-) = —ph 1 ^ are given by 

>!* = _ ( i + _L_ a a a ) (r» \Wfc . 
V n+l ) y 



-n 



P 



Note that, on the light cone (a — 0), p ~ ~{T v y ) and approach a constant value of order H n+2 . Thus, as expected, 
(Tjjj,) is completely regular on the horizon. 
For the + branch, we find 

n I 1 11 1 

irpy \bulk _ rrn+2 J ^ (o 07*1 

V y/ 25 (n+1) \{{aHf - 1 -nffr c {aH) 2 ^ 1 ) (1 - l/(aff) 2 )" k ' 

nHr c -2 2 (n+l-l/Hr c ) D /I 

Note that in spite of appearances, this expression is finite for i7r c = 1/n. At infinity (a ~ y — > oo), the leading 
order terms of the three summands in l|3.27|l cancel each other, so (T y y ) bulk and p buLk decay like H n+2 (ai/)~ 2 (™+ 2 ' ~ 
l/(Hy 2 ) ("+ 2 ). 

Close to the brane (y — + 0) , 

{n> w* K± ^L_*(i + , 

^<->(ra+l) y 

so (T v y ) hulk vanishes in the H — > limit, as required by the local conservation of the stress tensor. As for the 
components parallel to the brane (for n = 3), we find 

imi )b ulk 9 / 1 2 + 3ffr cT 9(ffr c ) 2 1 \ 

^ " 256^ 2 U 5 2(1 T ZHr c ) r c y* + ' ' ') ^ ^ * J 

The r c limit of this expression is not well defined because it holds only for y smaller than H and r c . 



C. Flat brane limit 



To recover flat brane case, we take £ ~ HAx where Ax is the distance between points x and x' , in n + 1 dimensional 
flat space. If we take i? — ► and keep -ffr c (and hence A>) fixed, this means that we are letting r c — > oo, and will 



We omit the anomaly term since it vanishes in the bulk for odd dimension, and on the brane it can be absorbed in a renormalization 
of the brane tension. 
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obtain the limit with Dirichlet boundary conditions. Recalling that in terms of the proper bulk coordinate y we have 
aH = e ±z = 1 ± Hy and taking the limit H — > in (|3~2UJl . lpT2T|) and 1|3.22|) . one easily obtains 

^ 1,(V ^ ;i ^ ) = -n5 ( ! + i,(l- + LT /2 " 

where Y = y + y' , which indeed corresponds to the renormalized Hadamard function in flat space with Dirichlet 
boundary conditions. 

To recover the H — ► limit keeping r c finite, we have to take into account that in this limit, 4 A2 — ► — 00. Using the 
representation 

1 f°° P - T 

-F\ (a, b, b, a+ 1; zi, z<x) = / dr 

a . n 



((1 - Zl e- T / a )(l - z 2 e- T / a )) b 



we obtain 



2 1 f e~ T 
lim -A G x = — / dr ^ . (3.29) 

°° nS {n+1) Jo (r 2 + (Y 2 + Ax 2 )/r 2 c - 2Y/r c ) n/2 



Specializing to n — 3, we obtain at Y = 0, 



y 77(1) , n n n 1 f 1 7T y 1 (Aa:/r c ) + g_ 1 (Aa:/r c , , 

limG (0,0;a:,a:) = — 5 --r- 3-3 r- , (3-30) 



H->0 



where Y„(z) is the Bessel function of the second kind and H v {z) is the Struve function. On the other hand, for Y ^ 
and Ax — 0, we have 

lim G (1) (y, y'; x, x) = — ^— — + — e Y ^T (l-n,-) 
yy,y ' ' nS (n+1) \ Y n r« V rj 



Note that in the limit r c — ► 0, we obtain G^ (y, y'; x, ar) = 1/ (niS^^Y™), recovering the result for an ordinary scalar 

with Ncuman boundary conditions. Conversely, in the r c — > 00 limit, we obtain G*' "* (y,y';x,x) = — l/(nSt n ±i\Y n ), 
corresponding to Dirichlet boundary conditions. 



As for the vacuum expectation of the stress tensor, one easily finds that (T v y ) bulk = as required by the local 
conservation of the stress tensor, and 



(T) 



i \bulk 



s, 



(n+l) 



(2y) 



n+2 



„n+2 



= 2y/r c 



r -l-n, 



2i) 



h' 



Near the brane, the stress tensor behaves like 



(T 



i \ bulk 
V 



256tt 2 



2 1 



2 1 



y 



r c y H 



+ 



4 1 



3 r 2 y 3 3 r\y 2 3 r^y 



8 log (y/r) 
3 r5 



(3.31) 



(3.32) 



The first term of this expression agrees with [23] and with the H 
beyond the first term holds for the self accelerated branch only. 
On the other hand, at large distances y r c , 



limit of Eq. Ij3.28|l . Note that the agreement 



(T 



i \bulk 
J I 



9 1 

— — o — k h\ (1 
256?r 2 y 5 1 y 



0(r/y)) 



(3.33) 



One observes that it behaves like close to the brane but with the sign flipped. This happens because the y ^S> r c 
{y -C r c ) regime is equivalent to the r c — ► (r — > 00) limit. In those limits, the field \& is subject to Neuman or 
Dirichlet boundary conditions respectively. As is well known, the sign of the leading (l/y 5 ) term depends on which 
condition is taken I28I. 



4 The contribution from Ai in I3.2UI vanishes in the limit H — » Q 
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D. Expectation value of the stress tensor on the brane 

The maximal symmetry of the y = const surfaces ensures that 

<a**a-*> = — !— (d k ^d k t>) hi . 

n + 1 J 

Hence, we can compute the v.e.v. of the surface stress tensor induced by the bulk held 4" as 

<T/) b ™" e = -(n - l)r c «%) = (n - l)r c ^ S(y) Urn ^—^ltlO_ . 

As is well known, the stress tensor on the brane diverges even after subtracting the Minkowski part. So, we have 
to introduce some alternative or additional method to regulate it. One option is to use dimensional regularization. 
Expressions like l/y n+2 are regularized by taking the y — ► limit for negative enough n, and then continuing back to 
n = 3. Using this prescription, we obtain 

ybrane _ 2 (n - 1) T (-1 - n) 4 _ _1_ ^ 1 1 log (r c fi) 

>(n+l) J-DT 



where we introduced an arbitrary renormalization scale /x. The divergence requires a renormalization of (p 4 , where 
<p = l/r c is the relevant held for the effective field theory, as argued in Section Hi Al 

To obtain a more physical picture, we shall show the form of the brane stress tensor regularized using a cutoff, 
which we can assume is related to the brane thickness d. To this end, we can evaluate the point split expression at 
y = and y' — d, and we obtain 5 



e d ^T(-4,d/r c ) 
2tt 2 " 

1/6 2 1 1 log(r c /x) 



(Ti) brane = ^ \r ,U/ ' CJ (3-35) 

S(y) h) (3.36) 



16tt 2 \d 4 d 3 r c d?r 2 c dr z c r 4 

where the dots denote positive powers of d and the expansion l|3.36|l holds for d <C r c . In the opposite limit, Eq. 
is perfectly regular, and in fact it vanishes for r c — > 0. Several comments are now in order. First of all, note that the 
finite parts of Eqs. I|3.34|l and (|3.36|l agree. Since Eq. I|3.36|l only holds for r c 3> d, we observe that the dimensional 
regularization result correctly captures the large r c form of (T fJiU ) brane , and should not be trusted for small r c . We 
shall interpret that the divergent terms in the expansion l|3.36|l should be absorbed in a (finite set of) counterterms 
involving negative powers of r c . This goes along the same lines of Sec. (|ll C|l . where we saw that the field associated 
to r c appearing in the 4D effective theory (in the large- r c regime) is nothing but tp = l/r c . 

Before discussing Eqns. H3.34[) - l|3.36[) anc ^ ^ e application of these results to the computation of the effective 
potential for r c , we shall display the form of the brane stress tensor when the brane is inflating. Again, dimensional 
regularization leads to a result valid for large r c . We shall not include here the form of (Tj) brane using a cutoff distance. 
The expressions are lengthy and do not illustrate essentially anything more than already contained in (|3.36|) . For the 
— branch, we obtain 

, , (n- l)H n+1 2 + nHr c ( 1 \ „, , s , , 

™° = S { ! +1) lT^Hf c B { n + 2+ H^>- n - 1 ) Siy)h > (3 - 3?) 
(1 + Hr c )(l + 2Hr c ) (2 + 3JTr c )(l + 4ffr c ) 



32tt 2 r 4 

+ ^ f 5 + + log (H/ri) - ^ + ...)> S(y) h l 

n-3 V Hr c ) &y /w (1 + 3Hr c )(2 + 3Hr c ) 1 



up to a rescaling of \x. 



5 Note that once we introduce a cutoff distance, in order to compute (T^ v ) we shouid use the full Green function instead of the subtracted 
one 13.161 . 
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For the + branch, we obtain 

, , (n - l)H n+1 2 - nHr c ( 1 . 

= (1 - Hr c ){\ - 2Hr c )(2 ~ 3ffr c )(l - 4ffr c ) 
32tt 2 r$ 



' V'f-E^--^ +log(ff/^)+ - — ^ Qtf ^ +...U(y)^ (3.38) 



n-3 r / (1 - 3Hr c )(2 - 3Hr c ) 

There are five types of divergent terms in l|3.37|l and l|3.38|l . proportional to H 4 , H 3 /r c , H 2 /r 2 , H/r 3 , and 1/rf,. 
They can be renormalized with counterterms of the type R%y K 3 ip, R^ip 2 , Kp 3 and ip 4 , where K is the extrinsic 
curvature and as before ip = l/r c . Note as well that we can readily identify two of the special cases when the divergent 
part of (|3.38|) vanishes: one corresponds to the truly self-accelerated solution (with r = 0, and Hr c = 1); the other 
corresponds to the only case when Eq. 1)2. 4|) has only one solution, Hr c — 1/2. 

In the H — ► limit, the finite part reduces to 



lD7T^ "" 4 



which agrees with the flat brane case for the self-accelerated branch only. This also happens in the expressions for 
the components of the bulk stress tensor Eqns. I|3.28|l and (|3.33|) . This might seem strange at first sight, but we shall 
this might be a consequence that in the — branch the bound state is present for all nonvanishing values of H, so this 
limit is discontinuous in this sense. 



IV. 1-LOOP EFFECTIVE POTENTIAL FOR r c 



This Section is divided into two parts. In the first, we obtain the 1-loop effective potential V e & for r c in the flat 
brane case, by three methods. The first one uses the expression that we obtained in Sec IIII Dl for the (expectation 
value of) the stress tensor on the brane. We re-derive it using the usual Kaluza-Klein decomposition and, finally, 
using the boundary effective theory of Section Til CI Then, we discuss the application of this potential to attempt to 
generate a large value of the cross-over scale r c in Sec IIV Al 



V c g obtained from the stress tensor on the brane, Eq. \3. 34t 

From (|2.8|l . the coupling between the (canonical) field tp associated to r c and \& is of the form r c (tp){d^) 2 . The 
equation of motion for ip is 

= r' c (p)(d*) 2 . (4.1) 

To take into account the quantum effects, one can substitute (S^) 2 by its vacuum expectation value. Then, we can 
identify the 1-loop effective potential V e g(ip c ) induced by W as 

V: s (r c ) = m) 2 ) , (4.2) 

and using Eq. H3.36J1 . 

VMr c ) = AA* + B^ + C^ + D^~ + . . . , (4.3) 

where A, B, C, D are constants, we identified A as the cutoff scale which we can associate to the brane thickness, and 
the dots indicate suppressed terms. As argued before, the relevant field in the 4D effective action is (p = l/r c . We 
discuss what renormalization conditions and the application of this result below. Before that, we shall re-derive the 
same result with two other methods. 
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Kaluza-Klein decomposition 

We can obtain slightly better understanding of the result i|4.3[) by working in terms of the dimensional or Kaluza- 
Klein (KK) reduction. Introducing the KK decomposition W — J dmU m 4> m as in Eq. (|3.2[1 . we see from Eq. 14. 2|) 
that 

Ks(rc)= I dmU 2 m ((diP m f) , (4.4) 



where ip m are 4D scalar fields with mass m. In the previous equation, we have assumed that modes with different 
m are independent, so that one has (ip m ip m ') = ('0m)^( m — m ')- We can easily compute ((ch/>m) 2 ) e.g. using point 
splitting. One finds 

{{d^mf) = ^ {16A 4 - 4m 2 A 2 - m 4 log (m/A) + . . .} , (4.5) 

where A is the cutoff and the dots stand for suppressed terms. For a Minkowski brane, the wave function of the KK 
modes on the brane is (see Eq. (|3.15[) ) 

1P\ = 1 1 

ml ° tt1 + (r c mf ' 

The integral over the KK continuum of (|4.5|l with this 'weight' contains a few divergent terms. The first term in 1)4. 5|1 
gives a quartic divergence, while the rest give at most cubic divergences. As we shall see below this is neatly captured 
by the boundary effective theory of Sec III CI Aside from this, we shall only comment that this results in a potential 
with the same features as (|4.3|l . In particular, we reproduce the same coefficient for the finite term. This agreement 
might not be surprising since, after all, the Green function (and so the stress tensor) have also been obtained as a mode 
sum. However, the regularization procedure is slightly different. Eq 1)4.3)1 is obtained by point-splitting the stress 
tensor, that is, evaluating the wave-function of the modes slightly off the brane. Instead, in the method described 
here, the sum over the modes is cut-off at a certain frequency. 



Boundary effective action description 

In the language of the boundary effective action of Sec 111 CI the decay of the resonant mode into higher dimensional 
modes is described by a (non-local) coupling of the form (py^Aip, where A = and <p — l/r c . Then, we can 

obtain the 1-loop effective potential as the sum over all 1-loop 1PI diagrams with vanishing the external momentum. 
One easily obtains 

1 f d 4 k f k/r c \ 



2 J (2tt) 4 ° V fc 2 

where k l is the Euclidean momentum flowing in the loop. This integral diverges. Introducing a cutoff A to regulate 
it, we obtain 

1 J4A 3 A 2 4 A log(r c A) + l/4 1 

^64^13-^ + R Vt + ' (47) 

where the dots denote terms suppressed by powers of A. 



A. A large cross-over scale? 

We shall regard Eq. 1)4. 3J1 as the result for the 1-loop effective potential V c s when the brane is flat (which holds 
for large r c ). The constants A,B, etc are unknown to the effective theory, and must be fixed by means of a set of 
renormalization conditions. We shall assume that ip — l/r c has a quartic potential aip 4 at tree level. Accordingly, 
we impose as renormalization conditions that the coefficients A, B, C and D in (|4.3(l vanish. Then, the potential at 
1-loop is of the Coleman- Weinberg [^J form 6 

V eS (<p) = cp 4 (a + /3 log^/A)) (4.8) 



We could also include the loops of ip itself, but the argument would not change much. 
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where (3 = l/64n 2 . The minimum of this potential is at r c ~ e a ^/A, which can be very large without fine tuning. 
This seems quite encouraging as it might provide a natural explanation why r c is so large in the DGP model, and 
ultimately why the current acceleration of the universe is so small. This potential gives a very small mass, of order 
l/r c . Aside from being potentially unstable under quantum corrections, this would require a very suppressed coupling 
of ip to matter in order not to generate an observable fifth force |36|. Rather than elaborating on these issues, we 
shall now address how this picture changes when we include the curvature of the brane - an essential ingredient if we 
are concerned about the accelerated expansion. 

For an inflating brane, the 1-loop contribution (which can be obtained from l|4.2|) and l|3.37fl . i|3.38fl ') is slightly 
more complicated than i|4.8[l . The most important difference is now the classical contribution ~ ±Hmpip discussed 
in Section Hi Al which clearly dominates over l|4.8|l . In the following, we shall assume that ip is not protected by any 
symmetry. It is then 'natural' to have a potential of the form 

Vesiip) = ±Hm 2 P ip + AA 4 + BA 3 tp + CA V + ■ • • , 

where A, B, ... are order-one coefficients and A is the cutoff. Given that the piece linear in ip represents the term due 
to the background curvature we shall take the renormalization condition B — (which amounts to a redefinition of 
mp). Then, the potential has a minimum (maximum) for the normal (self-accelerated) branch at r c ~ (A/mp) 2 C/H. 
There are different estimates available in the literature 0, ITol ITU IT^ . 13^ of the cutoff A for the bulk graviton, but all 
of them are much lower than mp. Hence, the expectation value generated for r c is too small for phenomcnological 
purposes. 

As an aside, we shall mention that an intriguing situation arises if the mass term for ip is Planckian. This is the 
case, for instance, if this field couples directly to matter and the cutoff in the matter loops is of order mp. One should 
still keep the linear term, as it encodes the coupling to the background curvature. Then, the extremum occurs at 
r c ~1/H, 7 which seems quite suggestive for phenomenology. This corresponds to a maximum for the self-accelerated 
branch and, in fact, it appears from Eq. 12. 4|) that generically there are no self-accelerated solutions with r c ~ 1/H. 
In the normal branch, instead, r c ~ 1/H is a minimum. This implies that effectively the value of r c 'tracks' the 
Hubble radius, modifying the whole setup. An analysis of the cosmology in this kind of scenario is more complicated 
(see [3(| for the case of a time-dependent Brans-Dicke field on the brane) and is left for future research. 

V. CONCLUSIONS 

We have discussed the quantum effects from the bulk fields in the DGP model, in a toy version of it where instead 
of the graviton we considered a bulk scalar field with a kinetic term on the brane. We have computed the Green 
function and the stress tensor induced by quantum fluctuations. We have shown that they are regular everywhere 
except on the brane, and described their form in several limits of interest. 

The cross-over scale r c in the DGP model is usually assigned a large value -of order of the present Hubble radius - in 
a somewhat ad hoc way. We used our results to clarify whether the quantum effects can provide a natural mechanism 
to generate such a large value. Once r c is promoted to a dynamical field, this problem translates to finding the 
potential that drives it to large values. Here, we considered the 1-loop effective potential V e s(r c ) induced by the bulk 
fields. 

With a flat brane in a flat bulk, we showed that this is possible by a quite simple argument. The fluctuations 
of the five dimensional 'graviton' provide a Coleman- Weinberg-type potential, which can induce a minimum at very 
large r c without fine tuning. When the curvature of the brane is included, the potential is largely modified due to a 
tree-level contribution originating from the background curvature. In the present setup, the value thus generated for 
r c is too small to be compatible with observations. However, simple extensions of the model suggest that it is possible 
to generate a minimum where r c tracks the Hubble radius 1 / H in the normal branch. 
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APPENDIX A: SOME FORMULAS FOR THE GREEN FUNCTION 

In order to perform the integral in (|3.17|l . it is convenient to use the representation for the massive Hadamard 
function in dS space 

(n- l)b( n ) (1 -cosC)t [ V 2 2 2 2 

T (% - ip)T (% +iv)T /l-cosCA.^ „ ( n n n + 1 1-cosC 

' 1 2-F1 -ip + iz,iV+ 77 



r(|_i P )r(i + ip)r(s=i) V 2 ) ^22' 2 ' 2 

where 2-Fi is the hypergeometric function and £ is the invariant distance between x and a;'. 
Including the bound state contribution, we obtain 

-1 / -1 \ h/2 oc oo \,\/ \ \ 



G q, ,r(n+l\ \ Ann') mil 1 2 



5 («) r (^) V 4aa '/ lEo^l A> -A< V"/2 + fc + j-A> n/2 + k + j-X 

(-l)T(« + 2fe + j) Kij2) TWH^) / l_cosC\ fc 

where f is the invariant distance in dS space and S7 n +i) = 27r 1+ " //2 /r(l+?i/2) is the volume of a unit n + 1 dimensional 
sphere. This sum can be readily performed for the first term in the curly brackets, leading to minus Go as defined in 
(|3.21|) . The two remaining terms in I|A1|) can be obtained from Go as 

G A = J da e Xa G (h -> e ±Q / 2 ff) , (A2) 

which leads to i|3.22[l quite immediately. 

For the sake of completeness, we include here the expression of the Green function when x and x' are 'aligned', 
that is, when their four dimensional coordinates coincide, £ = 0. We obtain 

g '" - - ^r i iw.i- + 2 ^ [a <a > ■ (wh2,ii) - a {x < ^' h ^ } ■ (A3 » 

where we introduced 

A(X,x) = Xx x B x g-A,l-n) 
and B x (a, b) is the incomplete Beta function. For small x, 

A{\x)^—^-x n ' 2 {l + 0{x)) . (A4) 
n/2 — A 

Let's now elaborate on the IR divergence present in (|A3|) for the — branch, due to the A> = n/2 contribution (the 
massless bound state). Using 

B x (a, b) = {-iy a B (1 - a - b, a) + (-l) b B 1/x (1 - a - b, b) 

where B(a, b) = B±(a, b) is the Eulcr Beta function, we can find that the expansion around A = n/2, 

nZ/2 1 

— — A (A, x) ~ - logx + (-l) n+1 B x (n,l — ri) - tp (n) - 1 + in + O (\ — n/2) , 

A A — n/2 
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up to an irrelevant constant, and ip(z) = V (z)/T(z) is the digamma function. In this equation, we kept the digmma 
term in order to obtain a finite limit for n = 3. For small x, 

(-l) n+1 B 1/x (n, 1 - n) - ip (n) ~ 7 - in + log a; + nx + ... 

up to constant terms and where the dots denote higher powers of x. The terms linear in logo; cancel, and the leading 
order term is O (x). The other two contributions to the Green function behave like const + O (x). Hence, the Green 
function is well behaved at infinity (which is represented by the light cone in the — branch). 
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